TWO TYPES OF PHASE SYNCHRONIZATION DESTRUCTION 
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In this article we report that there are two different types of destruction of the phase synchroniza- 
tion regime of chaotic oscillators depending on the parameter mismatch as well as in the case of the 
classical synchronization of periodic oscillators. When the parameter mismatch of the interacting 
chaotic oscillators is small enough the PS breaking takes place without the destruction of the phase 
coherence of chaotic attractors of oscillators. Alternatively, due to the large frequency detuning, the 
PS breaking is accomplished by loss of the phase coherence of the chaotic attractors. 
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Introduction 



Phase synchronization (PS) of chaotic oscillators ob- 
served in many systems (physical, technical, chemical, 
biological, physiological, etc.) has attracted great at- 
tention of scientists recently [H, 0, 01- PS has practical 
importance [J, |5|], with the phase synchronization con- 
cept involving theproblems of the phase coherence of 
chaotic attractor [y, 0] and occurrence (destruction) of 
PS regime 

The scientists usually do not distinguish between types 
of the phase synchronization regime destruction. As an 
exception we can mention work [9] where the existence 
of three types of transitions to PS regime in the cou- 
pled oscillators depending on the coherence properties of 
oscillations were described. Contrary to the conclusion 
obtained in [Q^], in this paper we report that there are 
two different types of the PS regime destruction depend- 
ing on the detuning of the system control parameters as 
well as in the case of the classical synchronization of the 
periodic oscillators. We find the following: that when the 
parameter mismatch of the interacting chaotic oscillators 
is small enough the PS breaking takes place without de- 
struction of the phase coherence of chaotic attractors of 
oscillators, whereas for the large parameter detuning with 
the decrease of the coupling strength between oscillators 
(or the amplitude of the external signal) the PS breaking 
is accompanied by the loss of the phase coherence of the 
chaotic attractors. To make our arguments more conve- 
nient we compare our methods and results with the ones 
described in work [9]. 

The structure of this paper is the following. Section [J 
presents two scenarios of the destruction of the phase 
synchronization regime in the unidirectionally coupled 
Rossler oscillators. Section [III explains the mechanisms 
resulting in these scenarios of the phase synchroniza- 
tion destruction. Section IIIII shows the relationship be- 
tween Lyapunov exponents and the phase synchroniza- 
tion boundary. Section [IVl displays the phase synchro- 
nization regime destruction in the systems with initially 
incoherent chaotic attractors. Section \V\ provides the fi- 
nal conclusion. 



I. DESTRUCTION OF PHASE 
SYNCHRONIZATION REGIME IN TWO 
UNIDIRECTIONALLY COUPLED ROSSLER 
SYSTEMS 

Let us start our consideration with two unidirectionally 
coupled Rossler systems 

Xd = -^dVd - Zd, 
yd = uJdXd + ayd, 

Zd =p^ Zd{Xd - c), 



Xq^ LUq^yrp Zf I £ (^X d '^7^)7 

yr = UJrXr + ayr, 

Zf P I Zf i^X rp C) , 



(1) 



where £ is a coupling parameter. The control parame- 
ter values have been selected by analogy with [iQ, [ll| as 
a = 0.15, p = 0.2, c = 10.0. The cj^-parameter which de- 
termines the main frequency of the response system has 
been selected as oor = 0.95, and the analogous parameter 
iJd of the drive system has been varied in the range from 
0.8 to 1.1 providing the mismatch of the interacting os- 
cillators. The mentioned above control parameter values 
provide the phase coherence of chaotic attractors of both 
drive and response systems for zero coupling strength. 

In Fig. [1] the phase synchronization area on the pa- 
rameter plane {cod'is) is shown. Following to Ref. Q we 
use two criteria to detect the presence of the phase syn- 
chronization regime. First, in the regime of the phase 
synchronization, locking of the mean frequencies 



^d 



d(pd{t) 
dt 



d(pr{t) 

dt 



(2) 



should take place, where (•) denotes the time average, ipd 
and (pr are instantaneous phases of drive and response 
oscillators, respectively. Second, the presence of PS may 
be detected by means of the examination of the phase 
locking condition [ll, 0, H, IH 



\^d{t) - (Pr{t)\ < const, 



(3) 



with the instantaneous phase of chaotic signal being 
introduced in a traditional way as the rotation angle 
(p = arctan(?//x). 
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FIG. 1: (Color online) The PS area on the (cj^, £)-plane. The 
PS boundary is shown in the range of the small parameter 
mismatch (the bold solid line) of coupled oscillators as well as 
in the range of the large one (the thin solid line). The chaotic 
attractor is incoherent in the region shown in grey color. The 
lower bound of the attractor non-coherence area is shown by 
the dashed line 



Examining the behavior of two unidirectionally cou- 
pled Rossler systems ([T]) we have found two different 
types of the PS regime destruction. When the parame- 
ter mismatch of the interacting chaotic oscillators is large 
enough, the PS breaking is accompanied by the loss of 
the phase coherence of the chaotic attractors, whereas 
for the small parameter detuning with the decrease of 
the coupling strength between oscillators PS breaking 
takes place without destruction of the phase coherence of 
chaotic attractors of oscillators. According to [1, 3] the 
attractor is phase coherent when the phase trajectory ro- 
tates around the origin of the projection plane without 
crossing it. Alternatively, the attractor is supposed to be 
incoherent. 

In Fig ^a,b the chaotic attractors of the response 
Rossler system and the dependence of the phase differ- 
ence (pd{t) — (Pr{t) on time t are shown below and above 
the onset of the PS regime for the small mismatch of co- 
parameters. One can easily see that with the decrease 
of the coupling strength e the PS regime is destroyed, 
with the chaotic attractor remaining phase coherent. For 
the large mistuning of the system control parameters the 
different scenario of the PS destruction takes place (see 
Fig. Elc^G?). Below the boundary of the PS area the 
chaotic attractor of the response Rossler system loses its 
phase coherence. 

Since the notion of the attractor coherence plays the 
key role in the chaotic synchronization theory the quanti- 
tative measure for the phase coherence should be defined. 
The measure of coherence may be characterized by means 
of the minimal distance p between the points of the phase 
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FIG. 2: (Color online) Attractors of the response Rossler 
system ([T]) and the dependencies of the phase differencies 
A(p{t) = ^d{t) — (pr{t) on time t for small detuning {cOd = 
0.91) of the drive Rossler system (a) above (s = 0.08) and (6) 
below {e = 0.075) the onset of the PS regime. Analogous fig- 
ures for the large {cod = 1.00) detuning of the cj^- parameter: 
(c) s is above {s = 0.127) and (d) below {s = 0.123) the PS 
boundary 



trajectory and the origin 



p = min ^/x^t)^^y^t), 



(4) 



where x(t), y{t) are the coordinates of the phase trajec- 
tory projection. Obviously, p > for the phase coherent 
chaotic attractor, and p becomes equal to zero when at- 
tractor loses its coherence. 

In Ref. [9] the coherence of chaotic attractor of the sys- 
tem under study was characterized by the sequence of the 
local maxima of y{t), i.e. max(?/). The chaotic attrac- 
tor is incoherent when some max(?/) < Yq occurs, where 
(Xo,lo) ^ is the fixed point around which the 

phase trajectory rotates. Such criterion is limited, since 
only one coordinate is considered. Our definition (j3|) of 
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FIG. 3: (Color online) The dependencies of the coherence 
measure p on the ^-parameter for the small (curve 1) and 
large (curve 2) system parameter mismatches. The arrows 
show the values of the coupling strength e corresponding to 
the PS destruction. 



the coherence measure allows to generalize the approach 
used in [9] in such a way for both coordinates to be taken 
into account. 

The dependencies of the coherence measure p of the 
chaotic attractor of the response system on the coupling 
strength e are shown in Fig. [3] for both small (curve 1) 
and large (curve 2) parameter mismatches. 

One can easily see that for small detuning of the control 
parameters (curve 1 in Fig. [3]) the measure of the coher- 
ence p is positive for all values of the coupling strength, 
whereas for the large parameter mistuning (curve 2 in 
Fig. [3]) there is an area just below the PS boundary where 
the coherence measure is equal to zero. Therefore, it 
is evident, that the chaotic attractor loses its coherence 
when the PS regime is destroyed in the case of large pa- 
rameter detuning and remains phase coherent both be- 
low and above the PS boundary if parameter mistuning 
is small enough. 

The area on the plane {uod^e) where the chaotic at- 
tractor of the response system is incoherent is shown in 
Fig. [1] in grey color. This region is joined closely to the 
boundary of the phase synchronization regime for the 
cjfi-values detuned essentially from the cjr-parameter of 
the response Rossler system. Therefore, we can come 
to conclusion that depending on mistuning of the con- 
trol parameters of system ([1]) the PS destruction occurs 
in two different ways, the mechanisms resulting to these 
types of PS destruction being supposed to be dissimilar. 

Therefore, the different manifestations of these mecha- 
nisms may be observed in the vicinity of the PS boundary. 
In particular, two different types of the pre- transitional 
behavior revealed below the PS boundary may be con- 
sidered as an indication of different mechanisms govern- 
ing the scenarios of the phase synchronization destruc- 



tion. Indeed, it was found that the type-I intermittency 
and the super-long laminar behavior (so-called "eyelet 
intermittency" [13]) take place for small differences in 
the natural freq uencies of the drive and response sys- 
tems [H, [3, [Hi m, 113 5 while, as far as large values of 
the natural frequency differences are concerned, the ring 
intermittency emerges [iS^. 

Let us show that the same both types of the synchro- 
nization destruction are also observed in the case of the 
driven periodic oscillators. This topic is discussed in de- 
tail in the next section. 



II. PHASE SYNCHRONIZATION 
DESTRUCTION OF VAN DER POL 
OSCILLATOR 

One of the approaches allowing to reveal different as- 
pects of the synchronization phenomenon is the consider- 
ation of the periodic oscillators (see., e.g. (1^J2Q|,[21L[22|). 
Therefore, let us use the classical model of the synchro- 
nization theory, namely, van der Pol oscillator 

X — {\ — x^)x + X = Asm.{(ji)et) (5) 

driven by the external harmonic signal with the ampli- 
tude A and frequency uOe to explain the mechanisms caus- 
ing the PS regime destruction in the coupled Rossler sys- 
tems ([1]). The value of the control parameter has been 
selected as A = 0.1. Oscillations in this case are certainly 
not chaotic but we can use the concepts of phase syn- 
chronization and phase coherence of an attractor in the 
same way as it is done for the chaotic oscillators. The 
phase of van der Pol oscillator is introduced as the ro- 
tation angle ip = arctan(?//x) as well as for the Rossler 
systems ([1]) considered above, whereas the phase of the 
external signal has been calculated as (pe{t) = ^et- 

The synchronization area on the plane {cOe ; A) is shown 
for the driven van der Pol oscillator ([5]) in Fig. [H As well 
as for the unidirectionally coupled Rossler systems ^ 
described in Sec. H two scenarios of the phase synchro- 
nization destruction may be observed. These two sce- 
narios are illustrated in Fig \5\ The attractors of van 
der Pol oscillator driven by the external harmonic sig- 
nal ([5]), and the dependencies of the phase difference 
A(f{t) = (f{t) — (pe{t) on time t are shown in Fig. [51 6 
for the small frequencies detuning (cjq — uJq) below and 
above the onset of the PS regime (here ujq = {^{t)) — 1.0 
is the natural frequency of the autonomous van der Pol 
oscillator). Evidently, with the decrease of the amplitude 
of the external signal A the PS regime is destroyed, with 
the attractor remaining phase coherent. Alternatively, 
for large mistuning of the frequencies, the different sce- 
nario of the PS destruction takes place (see Fig. El c^d), 
since the attractor loses its phase coherence below the 
boundary of the PS area. 

This statement may also be confirmed by means of the 
dependencies of the coherence measure p of the attractor 
of the driven van der Pol oscillator on the amplitude A 
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FIG. 4: (Color online) The PS area on the (cje, ^)-plane for 
the driven van der Pol oscillator. The PS boundary is shown 
in the range of the small frequency mismatch (the bold solid 
line) as well as in the range of the large one (the thin solid 
line). The attractor is incoherent in the narrow region shown 
in grey color. The lower bound of the attractor non-coherence 
area is shown by the dashed line 



of the external signal given in Fig. [6] for both small and 
large frequencies mismatches (cjq — ^e)- Again, as in the 
case of coupled Rossler systems (Fig. [3]) the measure of 
the coherence p is positive for all values of the amplitude 
of the external signal (curve 1 in Fig. [6|), if the frequen- 
cies ujQ and Ue are detuned slightly, i.e., the attractor of 
the driven van der Pol oscillator remains phase coher- 
ent both below and above the PS boundary. As far as 
the large frequency mistuning is considered (curve 2 in 
Fig. [6]) there is a narrow range located just below the PS 
boundary on which the coherence measure p is equal to 
zero, hence the attractor loses its phase coherence, when 
the PS regime is destroyed. 

The area on the plane (cje,^), where the attractor 
of the driven van der Pol oscillator ([5]) is incoherent, 
is shown in Fig. [H in grey color. As in the case of 
Rossler systems ^ considered in Sec. |I] this region is 
joined closely to the boundary of the phase synchroniza- 
tion regime (although the width of this area is less than 
in the case of Rossler oscillators) when the frequency of 
the external signal is detuned essentially from the natural 
frequency of van der Pol oscillator. 

Thus, when the PS regime is destroyed, the evolution 
of the attractor of van der Pol oscillator ([5|) driven by 
the external harmonic signal is the same as the evolution 
of the attractor of the response Rossler ^ system with 
a decrease of the coupling strength e, despite of the fact 
that the Rossler system dynamics is chaotic, and van 
der Pol oscillator one — periodic. Therefore, we can 
assume that the destruction of the phase synchronization 
regime is governed by the same mechanisms in both cases, 
with transition from synchronous to asynchronous state 
of van der Pol oscillator being possible to be explained 



-0.4 




-0.8 -0.4 0.4 X 



3000 6000 t 



FIG. 5: (Color online) Attractors of van der Pol oscillator 
driven by the external harmonic signal (|5]) and the dependen- 
cies of the phase difference A(p{t) = (p{t) — coet on time t for 
small (cje = 0.98) frequency detuning (a) above (A = 0.0250) 
and (b) below {A = 0.0230) the onset of the PS regime. Anal- 
ogous figures for the large {ooe = 0.9) detuning of the external 
signal frequency: (c) A is above {A = 0.0775) and (d) below 
(A=0.0750) the PS boundary 



analytically by means of the complex amplitude method. 

It is well known that there are two different scenarios 
for synchronization destruction in a periodic oscillator 
driven by an external force, corresponding to small and 
large detunings between the natural and external signal 
frequencies (see, e.g., tutorial [12]), respectively. Under 
certain conditions (i.e., for the periodically forced weakly 
nonlinear isochronous oscillator), the complex amplitude 
method may be used to find the solution describing the 
oscillator behavior in the form 



x{t) = Rea{t)e'' 



(6) 



For the complex amplitude a{t) one obtains the averaged 
(truncated) equations 

d = —iua -\- a — \a\'^a — ik^ (7) 
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FIG. 6: (Color online) The dependence of the coherence mea- 
sure p on the A-parameter for the small (curve 1, oje = 0.98) 
and large (curve 2, cje = 0.9) frequency mismatches. The 
arrows show the values of A-parameter corresponding to the 
PS destruction. In the frame the small part of the p{A) de- 
pendence near the point of the PS destruction for the large 
frequency detuning is shown 

where v is the frequency mismatch, and k is the (renor- 
malized) amplitude of the external force. For the small 
V and large k the stable fixed point on the plane of the 
complex amplitude 

a* = Ae^'^ = const (8) 

corresponds to the synchronous regime, with the synchro- 
nization destruction corresponding to the local saddle- 
node bifurcation associated with the global bifurcation of 
the limit cycle birth. Note, in this case the amplitude of 
the originated limit cycle is large initially. For the large 
frequency mismatches, the different mechanism of syn- 
chronization destruction is observed. With the decrease 
of /c-value the fixed point (stable node) on the complex 
amplitude plane becomes sequentially a stable focus and 
an unstable focus (via the Andronov-Hopf bifurcation) , 
with the limit cycle originating with the infinitesimal am- 
plitude. With the further decrease of /c-parameter, the 
amplitude of the limit cycle grows from zero monotoni- 
cally. In this case, the phase synchronization destruction 
is connected with the limit cycle location on the complex 
amplitude plane [12]. When the limit cycle starts en- 
veloping the origin, the synchronization regime begins to 
be destroyed. Obviously, in this case the attractor loses 
its phase coherence due to the fact that there are the 
moments of time when the modulus of amplitude \a{t)\ 
is equal to zero. As far as the small frequency mismatch 
is concerned, the attractor is always phase coherent due 
to \a{t)\ 7^ 0, Vt G (-oo; +oo). 

Since the use of Eq. ([6]) for application of the com- 
plex amplitude method is equivalent to the transition to 
the revolving coordinate system, two mechanisms of the 
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FIG. 7: (Color online) The phase trajectories of driven van 
der Pol oscillator on the rotating plane {x\y') for the small 
{a,b) and large {c,d) parameter mistuning, the A- values hav- 
ing been selected just above {a,c) and below {b,d) the PS 
boundary. Figures {a,c) correspond to the synchronization 
regime, figures {b,d) illustrate the asynchronous behavior 

phase synchronization destruction described above may 
be easily revealed by means of the consideration of the 
behavior of the driven oscillator on the plane {x' ^y') ro- 
tating with the frequency uOe of the external signal around 
the origin. These considerations of the rotating plane 
may be made apparent by using the coordinate transfor- 
mation 

x' = Xr cos (p -\-yr siu (f^ . . 

y' = —Xr sin (p -\- yr cos ^ ^ 

where cp = (pe{t) is the instantaneous phase of the exter- 
nal signal. 

In case of small frequency detuning {loq — uOe) the sta- 
ble node is observed on the rotating plane {x'^y') for 
the synchronized van der Pol oscillator (Fig. [7J a). As 
soon as the amplitude of the external signal starts to 
be below the phase synchronization boundary, the local 
saddle-node bifurcation, associated with the global bi- 
furcation of the limit cycle birth, takes place and on the 
(x', ?/')-plane the cycle with the initially large amplitude 
is becoming observable (Fig. [7J h) instead of the stable 
node. As far as the large frequency difference {uoq — cc^e) is 
concerned a cycle (may be smeared, if the A-parameter is 
not small enough) should be found on the rotating plane 
both below and above the PS synchronization bound- 
ary (Fig. [71 c^G?). With the decrease of the intensity of 
the external signal the amplitude of the cycle increases 
monotonically, with the destruction of the phase synchro- 
nization regime corresponding to the control parameter 
value when the cycle starts to cross the origin (Fig. [71 d). 
When the cycle does not envelop the origin (see Fig. [71 c), 
the phase synchronization regime takes place. 

Obviously, if the mechanisms causing the PS regime 
destruction are the same in the cases of both periodic and 
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FIG. 8: (Color online) The phase trajectories of the response 
Rossler system ([T]) on the rotating plane {x\ y') for small (a^6, 
ujd = 0.91) and large {c,d^ uoa = 1.0) parameter mistuning, the 
coupling strength values having been selected just above {a,c) 
and below {h,d) the PS boundary 



chaotic oscillators, one has to obtain the similar results 
for chaotic oscillators, too. Therefore, let us consider the 
dynamics of coupled Rossler systems ^ on the rotating 
plane ([9]) in the same way as it was done for the driven 
van der Pol oscillator, with the phase (f = (fdif) of the 
drive Rossler system being used instead of the phase of 
the harmonic signal. 

The results obtained for Rossler systems ^ are pre- 
sented in Fig. [SI For the small system parameter mis- 
match (Fig. [HI a) the behavior of the synchronized re- 
sponse oscillator in the vicinity of the PS boundary looks 
like noise smeared stable fixed point on the (x', 7/')-plane. 
This effect arises insofar as the Rossler system may be 
considered as a noise smeared periodic oscillator (see, 
e.g., [23]). When the coupling strength between Rossler 
oscillators decreases and the PS regime is destroyed, the 
stable noise smeared fixed point disappears and, in full 
agreement with the dynamics of the driven periodic oscil- 
lator (compare with Fig. [7J a,h)^ the noise smeared limit 
cycle with initially large amplitude occurs (Fig. ^h). 
Similarly, as well as in the case of van der Pol oscilla- 
tor, if the large parameter mismatches are concerned, the 
noise smeared cycle is observed both below and above 
the PS regime boundary (see Fig. [51 c^ti and compare it 
with Fig. [71 c,d)^ with the destruction of the phase syn- 
chronization regime corresponding to the cycle starting 
enveloping the origin (Fig. [H d). 

Thus, taken into account the results mentioned above 
we can draw a conclusion that the mechanisms of two 
types of the chaotic PS destruction are the same as in 
the case of the synchronization of periodic oscillations. 



III. LYAPUNOV EXPONENTS AND PHASE 
SYNCHRONIZATION BOUNDARY 

Let us discuss now the problem concerning the rela- 
tionship between the observed types of PS destruction 
and transitions to PS regime described in Ref. Q. In 
three different types of transitions from asynchronous 
dynamics to the phase synchronization regime were re- 
ported of, depending on the coherence properties of the 
chaotic attractors of the uncoupled systems. 

The degree of noncoherence of attractors of the au- 
tonomous oscillators was evaluated in [9] by means of 
the phase diffusion coefficient 



2D^t, 



(10) 



where (p{t) is the instantaneous phase of chaotic oscil- 
lator, (•) denoting the ensemble average. For a phase 
coherent chaotic attractor the phase increases approxi- 
mately uniformly and the value of D^^ is rather small, 
whereas for a funnel attractor the increase of the phase 
is strongly nonuniform and is a few orders larger in 
magnitude. 

As a criterion of distinction of different types of tran- 
sitions to/from the PS regime, the correlation between 
the dependency of the PS synchronization threshold (the 
critical curve li) and the dependencies of Lyapunov ex- 
ponents (LEs) (the critical curve I2 corresponds to the 
transition of one of the zero Lyapunov exponents to the 
negative values and the critical curve ^3 corresponds to 
passing through zero of the one of the positive Lyapunov 
exponents) on the coupling strength has been used. The 
types of transitions described in Q were the following. 

(I) Both uncoupled oscillators are characterized by the 
phase coherent chaotic attractors, therefore the value of 
the diffusion of phase is small (e.g., ~ 10~^ for 
the system considered in [9|). In this case, the phase 
synchronization regime is assumed to occur immediately 
after the transition of one of the zero Lyapunov exponents 
to a negative value, i.e., the critical curves h and I2 are 
close to each other, curve I2 preceding curve li. 

(II) Both chaotic attractors of the uncoupled oscillators 
are funnel, and the phase diffusion takes the intermedi- 
ate values in comparison with types I and III. In [9] the 



values of D,, 



IQ- 



- 10 were reported for this type of 



transitions. The critical curves li and I2 are supposed to 
be clearly separated, both lying below the critical curve 

(Ill) The uncoupled oscillators considered are charac- 
terized by highly incoherent chaotic attractors, and the 
phase diffusion is rather strong (for the system consid- 
ered in [9] Dy:, exceeds 10~^). In this case, the critical 
curve li is expected to lie above both curves I2 and ^3, 
i.e., the phase synchronization regime occurs after one 
of the positive Lyapunov exponents passes to negative 
values. 

So, Lyapunov exponents seem to be important to dis- 
tinguish between different types of transitions to/from 
the phase synchronization regime. Moreover, PS regime 
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FIG. 9: (Color online) The critical curves /i, h and I3 on 
the {(jJdjS) -plane. The critical curve h corresponds to the PS 
boundary, it is shown in the range of the small parameter 
mismatch (the bold solid line) of coupled oscillators as well 
as in the range of the large one (the thin solid line). The 
dotted line I2 corresponds to the transition of one of the zero 
Lyapunov exponents to negative values. The dashed line I3 
taken from our previous work [10] corresponds to zero crossing 
of the one of the positive Lyapunov exponents. The parameter 
plane is divided into three regions according to three types of 
transitions to phase synchronization described in [9] (regions 
I, II and III, respectively) 



arising is usually described in terms of Lyapunov expo- 
nents (LEs) (see, e.g. [ll, i, S, [o^ ) . Eventually, we have to 
compare our results concerning different types of transi- 
tions from/ to the PS regime with the ones given in 
Therefore, let us consider the relation between critical 
curves /i, I2 and Is corresponding to the onset of phase 
synchronization, the transition of one of the zero Lya- 
punov exponents to negative values and the passing of 
the one of the positive Lyapunov exponents through zero 
for the coupled Rossler oscillators ^ on the (ujd^s) pa- 
rameter plane, respectively (see Fig. [9]). 

As it was mentioned above, the autonomous sys- 
tems ([1]) under study with the pointed values of control 
parameters (and zero coupling strength) are character- 
ized by the phase coherent attractors in the whole range 
of cj-parameter variation. Correspondingly, the phase 
of chaotic signal increasing approximately uniformly, the 
phase diffusion coefficient D^^ is rather small and does not 
exceed the value of 10~^ for all values of cj-parameter. 
Therefore, we must expect that only the ffist type of 
transitions described in is observed for two unidirec- 
tionally coupled oscillators ([T]), whereas two other types 
(II and III, respectively) can not be observed for this 
system. Nevertheless, having considered the dynamics 
of two unidirectionally coupled Rossler systems ([T]), we 
have observed types II and III of transitions as well as 
type I (see Fig. [9]). 

As it was mentioned above, the transition of one of the 



zero Lyapunov exponents to negative values is supposed 
to be closely connected with the onset of PS in the case 
when coupled oscillators are characterized by the phase 
coherent attractors (see also [s, 9]). Our results refute 
this statement: Fig. [9] shows that the curve I2 coincides 
with the curve li (onset of PS) only in a very small range 
of the parameter mismatch of coupled oscillators (this 
range is shown in Fig. [9] in dark grey color and marked 
by number I) , although uncoupled both the drive and re- 
sponse Rossler systems are characterized by the phase co- 
herent attractors in the whole range of the a; (^-parameter 
variation. We can not say that the transition of one of the 
zero LEs to negative values determines PS arising even in 
the case of small mismatch of system parameters. Indeed, 
for the parameter cod = 0.9 separating two types of the 
PS destruction the coupling strength value Sps ^ 0.099 
corresponding to the onset of PS two times exceeds the 

-value corresponding to /2-curve {ei^ ^ 0.05) (see also 
Fig. [9]). Therefore, we come to conclusion that the tran- 
sition of one of the zero LEs to negative values precedes 
arising of the PS regime but does not determine it. We 
suppose that this transition may be connected with time 
scale synchronization (23, , but this problem requires 
further consideration. 

Thus, despite the fact, that both interacting chaotic 
oscillators are characterized by the phase coherent at- 
tractors for zero coupling strength the ffist type of tran- 
sitions described in Ref. Q is observed only in the narrow 
range of the oj^^-parameter values (see Fig. [9]). Moreover, 
two other types of transitions, described in [9} and sup- 
posed to correspond to the synchronization of oscillators 
with initially incoherent chaotic attractors, are also ob- 
served in the considered system ([1]) of oscillators with 
originally coherent chaotic attractors. Indeed, in Q it 
was reported that the case when the boundary of the PS 
regime arising (curve h) is between the critical curves I2 
and Is corresponds to the synchronization of oscillators 
with funnel chaotic attractors. Nevertheless, for the con- 
sidered system ^ this type of transitions takes place in 
the cc;^i-parameter value ranges shown in Fig. [9] in light 
grey color and labeled by number II. Similarly, the third 
type of transitions described in [9] and supposed to cor- 
respond to the systems with highly incoherent attractors 
(the critical curve li lies above both the critical curves I2 
and Is) is also observed in ([1]) despite of the initial coher- 
ence of system attractors (white regions III in Fig. [9]). 

It should also be noted that the location of the critical 
curve ^3 on the (cj^i; 6:)-plane coincides with the onset of 
the generalized synchronization (GS) regime [26l, [27I, [28| . 
The mechanisms determining the location of the GS 
boundary on the parameter plane have been considered 
for unidirectionally coupled Rossler systems ([1]) in our 
previous pape r llQj with the help of the modified sys- 
tem approach [11] . Based on the results given in [10] we 
can make a decision that the location of the GS bound- 
ary (and, correspondingly, the critical curve ^3) does not 
relate with phase synchronization phenomenon. 

Thus, we come to conclusion that for two unidirection- 
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ally coupled Rossler oscillators ([T]) two types of the PS 
destruction are possible in dependence on the system pa- 
rameter mismatch. Of course, one can distinguish three 
types of the PS destruction according to the location of 
the PS regime boundary and critical curves I2 and ^3 on 
the parameter plane, but there is no reason of doing that. 
Contrary to distinguishing PS destruction types between 
two classes (when there are well-known physical mecha- 
nisms determining such approach), the consideration of 
three types of transitions to PS based on the locus of 
the critical curves /i,2,3 on the parameter plane does not 
seem to be reasonable, since the positions of ^2,3 curves 
are caused by different mechanisms which do not concern 
the PS phenomenon. 



IV. PHASE SYNCHRONIZATION 
DESTRUCTION IN THE CHAOTIC SYSTEMS 
WITH INITIALLY INCOHERENT CHAOTIC 
ATTRACTORS 

Since the behavior of unidirectionally coupled Rossler 
systems having been examined we came to conclusion 
that for the oscillators with initially phase coherent at- 
tractors ([1]) two types (rather than three) of the PS de- 
struction should be distinguished, it is quite reasonable 
to learn what types of PS destruction may be observed 
if two oscillators with initially incoherent chaotic attrac- 
tors are considered. To clarify this point we have studied 
two mutually coupled Rossler systems 

^1,2 = -^1,2^1,2 - ^1,2, 

yi,2 = ^1,2^1,2 + tt^l,2 + d{y2,l - ^1,2), (11) 

i^i,2 = 0.1 + 2:1,2(^1,2 - 8.5) 

from the point of view of the concept discussed in Sec. [J 
andini Note, it is the system that were considered in 
for which three types of transitions connected with the 
attractor coherence properties and critical curves loca- 
tion had been described (Fig. [TQ|) . 

The control parameter values of the systems (fTTI) have 
been selected the same as they were given in [9]. In 
Eq. ([TT]) d is the coupling strength, uoi = 0.98, 002 = 1.02, 
the parameter a G [0.15; 0.3] determines the topology of 
chaotic attractors. It is known, that when a exceeds the 
critical value ac {ad ~ 0.186 for coi and ac2 ^ 0.195 for 
UJ2) the chaotic attractor of autonomous Rossler system 
becomes incoherent [9j. Following [9.] the phase ip{t) has 
been defined as the rotation angle (p = arctan(7//i:) on 
the plane (x; y) to study PS of the systems with initially 
funnel attractors. 

Having studied the behavior of two mutually coupled 
Rossler systems (pT]) along the boundary of PS regime 
on the plane (a^d) given in [9], we have also found 
two scenarios of the PS regime destruction described 
in Sec. [J as well as for the case of two unidirectionally 
coupled oscillators ([1]). For small values of a-parameter 
(a < a* « 0.205) the first type of the PS regime destruc- 
tion takes place when chaotic attractors of both systems 




0.00 ' ' — ' ' ' ^ 

0.15 0.20 0.25 0.30 

FIG. 10: The phase synchronization boundary of two coupled 
Rossler systems (pT]) . The bold solid line corresponds to the 
phase synchronization regime destruction when both chaotic 
attractors retain their phase coherence both above and be- 
low the synchronization threshold. The thin line presents the 
synchronization boundary when the PS destruction is accom- 
panied by break of the attractor coherence. The change of the 
type of the PS regime destruction takes place at a* ~ 0.205. 
The attractor of the first Rossler system is incoherent in the 
region shown in grey color. The lower bound of the attractor 
non-coherence area is shown by the dashed line. The areas 
corresponding to three types of transitions to PS regime de- 
scribed in [9] are labeled by numbers I, II and III, respectively 



keep their coherence on the plane (x, y) both above and 
below the boundary of PS (Fig. [TO]). When a > the 
chaotic attractor of the first system loses its coherence 
as soon as the coupling strength d is below the onset of 
PS, with the attractor of the second system remaining to 
be phase coherent. Note, that a* > ad, 2 therefore the 
change of the PS destruction type takes place when both 
coupled systems are characterized by the initially funnel 
chaotic attractors on the plane (x^y). 

To make our decision more convincing, we have also 
calculated the coherence measure p vs. coupling pa- 
rameter d for two a-parameter values corresponding to 
both types of the PS regime destruction, to be exact 
a/ = 0.19 < and a// = 0.21 > a* (Fig. [TT]) . Since the 
evolution of the chaotic attractors is considered on the 
velocity plane (x, y) we have used 



p= min ^/i^t)Tfit) (12) 

instead of Eq. (|4]) . Evidently, the chaotic attractor of the 
first Rossler system remains phase coherent both below 
and above the PS boundary for a/ = 0.19 and loses its 
coherence when the PS regime is destroyed in the case of 
an 0.21. 

The evolution of the chaotic attractor of the first 
Rossler system, when the PS regime is destroyed, is 
shown in Fig [121 The first type of the PS regime destruc- 
tion, when the chaotic attractors are phase coherent both 
above and below the PS boundary, is shown in Fig[T2l a,b. 
Fig- [El c,d illustrates the second type of transitions ac- 
companied by the destruction of the coherence of one of 
the chaotic attractors just below the PS boundary. 

Thus, for two coupled systems with initially funnel at- 
tractors, the phase coherence of chaotic attractor on the 
velocity plane (x^y) is destroyed in the same way as it 
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FIG. 11: (Color online) The dependencies of the coherence 
measure p on the coupling strength d for ai = 0.19 (curve 1) 
and an = 0.21 (curve 2). The arrows show the threshold of 
the PS regime 



was observed on the plane (x, y) for the systems ([T]) with 
initially phase coherent attr actors. 

Taking into account all obtained results we come to 
conclusion that for the system of two mutually coupled 
Rossler oscillators with initially phase incoherent attrac- 
tors (pT]) as well as for the oscillators with phase coherent 
attractors ([1]) two types (instead of three) of the PS de- 
struction should be distinguished. This conclusion is also 
confirmed by Fig. 6 given in [9] where the mean frequency 
ratio Vti/Vt2 versus coupling strength is shown. The be- 
havior of Vti/Vt2 is quite different for the first and the 
second types of transitions to PS proposed in [9] but is 
the same for the second and the third ones. 
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FIG. 12: (Color online) Attractors of the first Rossler sys- 
tem (pT]) for two values of a-parameter corresponding to two 
types of the PS regime destruction. Figures {a,h) show the 
evolution of the system behavior on the velocity plane (x,y) 
for ai — 0.19 when the attractor remains phase coherent both 
above (a) and below (6) the PS regime boundary. The second 
scenario of the PS regime destruction {an = 0.21) accompa- 
nied by the phase coherence loss of the chaotic attractor is 
shown in Fig. [c,d). The behavior of coupled oscillators is 
synchronized in Fig. (c) when the attractor of the first Rossler 
system on the plane {x^y) is phase coherent. Alternatively, 
the phase incoherent attractor is shown in Fig. (c?) after the 
PS regime has been destroyed. In the frames the trajectories 
on the plane {x^y) near the origin are shown. 



V. CONCLUSION 

We have described two different types of the PS regime 
destruction as well as the mechanisms resulting in the 
destruction, which are the same as in the case of the 
classical synchronization of periodic oscillators. The first 
of them is caused by the loss of the common rhythm 
of chaotic oscillations and the second one is caused by 
the loss of phase coherence of chaotic attractor. These 
types have been observed in the systems with initially 
both phase coherent and funnel chaotic attractors. The 
transition of one of the zero Lyapunov exponents to neg- 
ative values precedes arising of the PS regime but does 
not cause it is the other important result. Similarly, the 
passing of the one of the positive Lyapunov exponent 
through zero does not relate to PS phenomenon. There- 
fore, there is no reason to use these LEs for the PS onset 
description. 

We assume that the number of the possible transi- 
tions from the phase synchronization regime to the asyn- 
chronous dynamics might not be limited by two types 
only which are described in our paper. For the more com- 
plex systems (e.g., for the system with the larger dimen- 



sion of the phase space) the breakdown of PS may also be 
associated with global bifurcations, e.g., homoclinic or- 
bits, Takens-Bogdanov bifurcation points, etc. [29]. Nev- 
ertheless, we suppose that the two considered types of 
the PS regime destruction should be typical for a wide 
class of nonlinear systems, e.g., such as Pierce diode [30| 
or laser system [16]. In particular, we have observed the 
same types of behavior in the vicinity of the PS bound- 
ary for the case of two unidirectionally coupled generators 
with tunnel diodes described in Ref. [31]. 
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